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The study of linear perturbation tlieory for general functions of the Ricci and Gauss-Bonnet 
scalars is done over an empty anisotropic universe, i.e. the Kasner-type background, in order to 
show that an anisotropic background in general has ghost degrees of freedom, which are absent on 
Friedmann-Lemaitre- Robertson- Walker (FLRW) backgrounds. The study of the scalar perturbation 
reveals that on this background the number of independent propagating degrees of freedom is four 
and reduces to three on FLRW backgrounds, as one mode becomes highly massive to decouple from 
C ' the physical spectrum. When this mode remains physical, there is inevitably a ghost mode. 

00 ■ I. INTRODUCTION 

<N . 

Modifications of gravity have been used in many contexts [l| . Perhaps the most famous is the first model of inflation 
, introduced by Starobinsky in 1980 01- However, in the last few years, modifications of gravity have been attracting 
' attention as an alternative to quintessence Q, in order to explain the late-time acceleration at large cosmological 
• , scales 0-01 • The most general modifications of gravity whose Lagrangian is built using only the metric tensors so as 
' not to introduce any extra vector or spin-2 degrees of freedom other than the graviton are such that the Lagrangian 
is given by 
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where R and G = — 4:R^i,R^" -'rR^uapR^'^'^^ are the Ricci [R) and Gauss-Bonnet [Q) scalars, respectively. We have 
fSj ' only two non-zero Lovelock scalars in four dimensions Q ■ There are several recent studies on this type of models 0, 0] . 
^ ' In Ref Q the linear cosmological perturbation theory of these general modifications on a Friedmann-Lemaitre- 
,. Robertson- Walker (FLRW) background was studied, and the result is that only one propagating mode for the scalar- 
type perturbations is present for a general function /. Furthermore, such a mode has a scale-dependent speed of 
propagation, depending on the wave number k as well as the time-dependence through the background quantities. 
Furthermore, the background can be chosen so that this scalar mode is not a ghost. Each term in the Lagrangian 
of this scalar mode has only two temporal differentiations, whereas the square of the Laplacian operator appears, 
leading to a dispersion relation such as = Bk^/a^, where i3 is a function of the background quantities. Even 
though the Lagrangian looks exotic, there is no ghost as long as the signature of the second temporal derivative term 
is normal [Toj . 

However, a few questions naturally arise by considering this result carefully. For example, in Q the studied 
background was special, i.e. spatially homogeneous and isotropic. If we consider other inhomogeneous/anisotropic 
\^ • backgrounds, the term proportional to fc*, i.e. the operator A^, may shift to the term proportional to uj^, i.e. the 
9^ ! operator df, which would imply the existence of ghost degrees of freedom. Therefore it is interesting to study the 
behavior of this same theory on one of the simplest anisotropic backgrounds, the Kasner-type solution. 
There is another point which is tightly connected to the previous one. The action ([1]) can be rewritten as 

r 

^ I d^xV^[FR + ^g-V{F,0], (2) 
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where F, and ^ are auxiliary fields. By integrating out these auxiliary fields, we can verify that the action ([2]) 
is equivalent to ([T]) [ll|. On the FLRW background among the two extra scalar fields introduced in Eq. ([2]), one 
linear combination of their perturbations does not have kinetic term. Here, it is totally unclear if such a feature is 
common to all backgrounds or it only happens on the FLRW one. In fact, as we shall see later, both modes propagate 
independently on Kasner-type backgrounds. In this paper, we will also see that having one more propagating field 
implies that the cx fc* dispersion relation returns back to a normal one, i.e. oc fc^. But how can we explain this 
apparent reduction of the propagating degrees of freedom from one background to another? As we shall see later, 
the determinant of the kinetic operator for the perturbation modes reduces to zero as we take the limit of the FLRW 
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backgrounds. At the same time, the mass for the mode whose kinetic term vanishes on the FLRW backgrounds in 
general blows up to infinity. Therefore, in the FLRW limit this extra mode decouples from the physical spectrum to 
be integrated out, and the number of effective propagating modes is reduced by one. 



II. BACKGROUND MODELS 
A. Kasner-type backgrounds 

For the reasons mentioned above, we study a simplified Kasner-type anisotropic background, which can be written 

as 

ds2 = -dt^ + /i (t) dx^ + /2 it) (dy2 + dz2) . (3) 

We analyze the behavior of the perturbations at linear order, by using 1 + 1 + 2 decomposition according to the 
symmetry of the background metric. From the symmetry in y—z plane, in general, the perturbation can be decomposed 
into three even-parity modes and one odd-parity mode. One of the three even-parity modes has a vanishing kinetic 
term in the FLRW limit, where the mode decouples since its mass becomes infinitely large. 
The equations of motion for the system ([2]) are 

dV dV ^ 

F G^u ~ , (5) 

where E^^ is the tensor defined by 

V^V.F - .g^.DF + 2i?V^V,e - 2g^.i?ne - 4i?^^VAV,e 

-4i?/VAV^e + 4i?^.ne + 4g^.i?"'^V„V^e + 4i?^a^.V"V'^e - ^ 9^..V . (6) 
On the Kasner-type background we have 

R 2 /2/1/1/2 + 4 hfif2 + 2 A/i/i - fifi ~ fifl 

mr? ' 

/2 f 2 /2/1/2/1 - 2 /i/lA + 4 /2/1 /1/2 - /2/2A' 



(8) 



We can use the (tt)-component of Einstein equations in order to write V{t) in terms of the other background quantities 

as 

2 f /2/2 A + 6 e/l A + Ffih + 4 f /2/2/1 + 2 fA/I 

^^'^ Wfi ■ 

Furthermore, one can use the (xx)- and (?/?/)-components of Einstein equations in order to write F and ^ in terms of 
the other background quantities as follows: 



^ = ,2,i\ 7TT f2flf2hF + 2 flflhP - 2 hfiFhfl - 4 fUhf! + 4 flMfi - hflf^F 

^hfl (/2/1 - /2/1) 

-fIfMF + 4 A/ie/2A A - 4 hf2ififi - 2 Ffffi + 2 AVi'^^ + 4 hfifihF] , (10) 

^' = "TTTTTT^T 7TT FhfH! + 2 i^A./IA - 2 /i^A/l + 2 A^A/| + 2 A^AA'^ + A^/l 

4/2/1/2 (/2/1 - /2/1) 

+8 .AeAA/2 - 4 /fe/l - 2 /ii^A./l - Fh.fi hh - 4 /i./|^A A - 4 /i/leA A] • (11) 



B. Almost FLRW backgrounds 



For practical construction of cosmological models, we are particularly interested in the behavior of these theories 
on the backgrounds close to FLRW. Therefore it is sometimes convenient to consider Kasner-type backgrounds which 
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are described by homogeneously perturbed FLRW as 

h=a{tf + 5Mt), (12) 

f2^aitf+Shit), (13) 

F^Fo{t) + SF{t), (14) 

^^m+m)- (15) 

At the hnear order in this perturbation, eombining the (xx)- and (?/?/)-components of the Einstein equations leads to 
a elosed differential equation for the field 

6 = 6h- 5f2 , (16) 

whieh reads 

S + 2i^S + M^S = 0, (17) 



where 



>^it) = -lH+ ^ ±(Fo + mo). (18) 

2 2{Fo + 4H(o) ™ 

M\t) = -2H-2H^- A(f„ + 4geo), (19) 

with H :— a I a. The differential equation admits the general solution 

b = cxa{tf \ : +C2a{tf. (20) 

' ^ W a(t')MFo(i') + 4i?Co(i')] 

The second solution represents just the constant shift of the scale of spatial coordinates. Namely it is just a gauge 
mode. The first physical solution decays, compared to the term, when the term a^(Fo + 4i7^o) increases in time 
faster than . In this case /i gets closer to /2, and the Kasner-type background reduces to FLRW. 

III. PERTURBATION ANALYSIS 

A. Reduction of the action for the even-parity modes 

We begin with the even-parity modes, which here mean the scalar-type perturbations in the two-dimensional space 
spun by y — z plane. In order to achieve this goal we consider the perturbation given by [l^l 

F ^F(t) + 8F{t,x), (21) 

i^m^m.^). (22) 

ds^ = -[1 + a(t,a;)]dt2 +a^;;^(t,a;)dtd.T-|-5j,Cdidy-(-a3Cdidz-^ [/i -H/3(t,a;)]da-2 + /2(t) (dy^ + dz^) . (23) 

We have fixed the gauge so that the perturbations of the spatial metric in the two dimensional subspace spun by y 
and z vanish. We will find it convenient to make a field redefinition. 



V-^JF+^^e, (24) 
il 



Expanding the action to the second order, and introducing a Fourier expansion for the fields with the wavenumber 
fc, we eliminate the perturbations of lapse and shift (a, x ^-^d C) by using the constraint equations obtained by the 
variations of the action with respect to these non-dynamical variables. Then, one reaches the result 



dtd3fc(/l/|)l/2 - B,j<\>,^j - C^,(j)^(bJ - D,j(t)^(j)j . (25) 



where (pi = {j3,ip,5^) and we have introduced q = {ky + k^Y^^. All the matrices A,B,C and D depend on the 
momentum only through the ratio 7 = k^/q. Only the matrix B is antisymmetric, while the others are all symmetric. 
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B. Degrees of freedom and the signature of kinetic term 

It is important to establish the sign of the eigenvalues of the kinetic matrix A. If they are not all positive, it means 
that there exist ghost-like degrees of freedom. The first key issue is that the matrix A is not degenerate on general 
backgrounds. Namely, 

does not vanish. Here we defined 

Ai = /i^^ + 2/i^, (27) 

J = 2Fj^fl + 2Fhh + fiFf2 + Gl^fU + 2Fh^^f2 + 6/1C/2 + A^^/i ■ (28) 

Thus, the number of independent degrees of freedom present in the even-parity perturbation is three. 

We notice that det A vanishes in the FLRW limit, /i — >■ /2- This implies a reduction of the number of the degrees of 
freedom in this limit. In fact, in this limit both ^13, and A23 become of order 0(6); whereas ^33, like the determinant, 
is of order 0{S'^), as shown by the following expressions 

Ai3 = : : : (~AOH^F£, -AHFp-f^ - 48i/V^^ - 144 TJ^^ 

a5 (F + 4 i^O (72 + 1) (2 Fi? + F + 12 ^iJ2)2 ^ 

-8 n^F-i^i - 6 HFF + ^2^2 _ 2 H^F'^ - F^ - 16 H^F^i^i - 32 FH^Cl , (29) 
^ _ 4ci (-8 Fj^H - 24 jj^H^ + 2 j^F ~ 36 jR'^ -8FH -F) 
~ {2FH + F + 12 ^7?2)2(^2 + 1) + 4 ^3 ' 

1 g 2 rr2 

A33 ^ : i : : (F^ ^ _ //S^ 2 t _^ 4 HFFj^ - 16 FH^-f^f 

a^{F + iH^f{-/^ + iy{2FH + F+12^H^)^ 
+2i H^F-f'^i ~ 48 H'^j'^i^ - 4 H^-f^F^ - HFF - 4 FH^i - 8 F^j^H^ - 44 H^F^ - 5 H^F^ - 96 H'^i^ 
+4 HFF-f'^ + 16 H^F-f^O , (31) 
(recall that ci is the magnitude of 5 = Sfi — 8/2), whereas the determinant tends to 

det A = : i : : . 32) 

aio(F + 4i/0(72 + l)2(2Fi? -f F 12^i?2)2 

The fact that the sign of dot A is opposite to the one of F + 4iJ^ shows that on Kasner-type backgrounds close to 
FLRW models, at least one mode becomes a ghost as F + 4iJ^ > is required for the tensor modes not to be ghost 
in FLRW backgrounds 

The matrix A, in the FLRW limit (ci 0), becomes block diagonal, i.e. it is composed of a 2x2 matrix and 1x1 
matrix. The latter is zero valued. The 2x2 submatrix is non-degenerate and finite. The sign of its eigenvalues are 
read from 

A22 = ^^^^X±^, (33) 
(2Fi/ + F+12ei?2)2 ' 

^11^22 - AL = ^ : . (34) 

(72 + l)2(2Fi/ + F + 12^i/2)2 a4 

Since both of these two factors are positive (assuming F -f AH^ > 0) , the (unique) ghost mode in this limit is identified 
with the third field, which is (5^ in the present calculation. We note that the important point to obtain the property 
that the kinetic term of the third field vanishes as is in how we choose the first and the second fields f3 and ip. As 
we will see below, the ghost mode S£, can be integrated out from the action near the FLRW backgrounds. 

C. The Laplacian operator and the behavior on FLRW 

Let us restrict our attention to better understand the FLRW limit for the Laplacian matrix C. In particular we 
are interested in the elements regarding the coupling of SS, with itself and with the other two fields. We find 

,,,C^-'-^^^^±m±J^l^, (35) 
(2FH + F + 12^2^)2 flio 
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and 

C,3 = + . (36) 

{12^H^ + 2FH + F)a'^ 

{12im + 2FH + F)a^ 
C 64ciij(27^ + 1)H^ 

{12im + 2F H + F){F + AH S,)a^' ' 

Therefore, also the g^-self-couphng vanishes on FLRW, but the g^-couphngs with the other fields do not. It is of 
crucial importance whether the remaining term quadratic in S^, the mass term vanishes or not. D^t, is given by 

^ ^ {F + 4^HY d^V 6{F + 4^H^) ^ 48HH^ 768H^H^{F + AHQ 

which does not vanish in generaQ- This implies that in the Lagrangian in the FLRW limit the only term quadratic in 
(5^ is its mass. It should be noticed that in constructing cosmological models we have to require _F w 1 at early times. 
Then, we have F H and <C 1. This implies that D^^/H'^ in general becomes very large at early times. This 
fact may help integrating out the ghost, as it naturally gives a very large mass to compared to its kinetic term. 
The part of the Lagrangian related to is written as 

C3-^a^D33 - a^S^ {q^ C13 /3 + Bi3/3 + C23 + B23V') + • • ■ + 0{5) . (40) 

Integratiiig out (5^, with the aid of the relation C13B23 = C23i?i3, we obtain the terms proportional to q^, which were 
found in tsl. 



D. Odd-parity modes 



Let us consider the odd-parity modes, i.e. the modes that are composed of divergence-less vectors in terms of 
two dimensions spun by y and z. Perturbations of all four-dimensional scalars are zero for odd-parity modes. The 
perturbation of odd-parity modes appears only in the metric perturbation as 







(41) 



where i,j = y or z, and Vi and Wi are divergence-less vector. We used one gauge degree of freedom belonging to the 
odd-parity modes to eliminate the perturbation in (ij)-components. 

As in the case of even-parity modes, we erase the shift perturbation Vi by using the constraint equation obtained 
by taking the variation with respect to Vi itself. Then, we obtain 



dM3fc/yv2A(„ 



Ji 



(42) 



where 



A, 



A1A2 



/i/|(Ai+72A2) ' 
(l+72)-i(Ai-f7'A2) 

A2Ai/2(/2/l-/2/l)/2 



{~2Fhfif! + Fflfih + 2F/1/1/I - flFfl - 2Fhflfl 



-ifUfI ~ ^fUhf2f2 + ififUhfi - 2/1^/1/2/2 + ^fifUfih + ^fiFfifi - 2fUf2fi) 



(43) 



(44) 



m 



i-Sf^FhftMfi + Sf2fiftef2ft - ^Fhf^ftYU2 - ^f^FfiftYUi + ^F'ftYfif2fi 



^ The expression ^ gpi ^ '^^^ rewritten in terms of second derivatives of / in action (1), as follows (^ gjX ^ ~ fss/ifRRfsS ~ fiig)' 



where fRR = (l^)gg 
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-8f2flfh'^'fi ~ ^khlHflFfl + SfififWifiF + ^FflhfWhF - GFhfifh'ih 
+12/1^/2/2^/2/17' - ^FflMii'Mfi - 2Fh-i^flFfl + 2Fflf2iHfl - 2Fflflf!j'F 
-^f2flfWiF - iF^fififWh + 8hfU'hf! + 2F^flflhh - 2F^flhflh + F^fh^fl) 
x[iFh + 2M)iFh + 2i-/^f2 + 2M + /2F7')(/2/i - /2/i)/i/|]"' , (45) 

and A2 = /2 + 2/2 ^. Here, C(-o) looks divergent in the FLRW limit, but in fact it is not. In the FLRW limit we find 

2 f + 4| 

2 2(8HHi + HF + 4:H^C + H^F + HF + 4:H^h , , 

m/_-, = ^ : . (48) 

w F + AH^ 

and Cy reduces here to the speed of the tensor modes in the FLRW background. 



IV. IMPLICATION OF THE PERTURBATION ANALYSIS AROUND KASNER-TYPE 

BACKGROUNDS 



A. Inevitable ghost 



Here wc show that the no-ghost condition > with ^(o) given in Eq. ()43p is incompatible with another no-ghost 
condition detA > for the even-parity modes. From Eq. (|26p . we find that Ai < is required for the absence of 
ghost. Then, if A2 < 0, A(^o) becomes negative for any 7, and hence the odd-parity modes becomes ghost. Even 
if A2 > 0, becomes negative for a sufficiently large 7, and hence at least a part of odd-parity modes becomes 
ghost. Therefore we conclude that the presence of ghost is inevitable for generic models of f{R, G) once we consider 
Kasner-type backgrounds. This is the main claim of this paper. 

However, as we have seen above, the perturbative action for such a mode does not completely vanish in general 
even when the kinetic term vanishes, e.g. the matrix D is not degenerate even in the FLRW limit. Therefore this 
ghost mode does not becomes strongly coupled in the FLRW limit, but it becomes infinitely massive. Thus, such a 
mode is to be integrated out as long as backgrounds are chosen to be sufficiently close to the FLRW universe. In this 
sense, such models might be harmless in spite of the presence of ghost. 



B. Particular cases with only three degrees of freedom 
We will consider here the case for which one has the condition 







diF,0 




\dRdgJ 





(49) 



In this case perturbations of F and ^ are not independent. For generic backgrounds, this is the case if £ = F{4))R + 

Below, let us consider special cases in which /(i?, t/) is a function in the form oi f{R + Xg) with A being a constant 
In this case, 



5^ 



(50) 



generally holds. Therefore one can eliminate 5£, from the perturbation action so that the number of even-parity degrees 
of freedom is reduced to two. Then, the no-ghost conditions for even-parity modes reduce to only two conditions 

^'^t ^ = - %tttt (2/iC + hF){Ffj + e/l) 
F^ftf2J^ 
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X i2fiFM - ifiFFh - 8/2/2^/1 - 4/2/2ei^/i - Qf^Ue + /2 ^^/iC) > 0, (51) 
fiF^fiJ^{Ff2 + 2^f2) 

~6f2M'Ff^Ffi ~ fififfP'F' - UfifU^FhFfl - ifihiPflF^fl ~ iflfiF'f^Fifi ~ flUP^hF' 
-f2hf2f?F^F' - Sf^f^F^M'fi - QfVWFFf^ - 2f2fU'fiF'f^ - 16f2fU'Fhf^ + 2eVi7'^/J 
-^ftflF^fiFik - 8f2^'Fflf! - AfU'fiFfl - SfifU'Ffi - Wfifii'F^'f^ - UhM^'l^Ffl 
-8fiM'F-f\fl - 22f2fUWFfl + f2i^hF^fh^ - 2flflflF^Fi - AflM^F^f^fi - mf2i'FfJlj^ 
-^flF^hl'F'fl - if^fiflF'j'F' - IGflF'hl'^'f^ - niFF^hfl - ^flfU'FVh' " ^flM'F'fh' 
+8f^f?F^^^^fi + iflfiF'j^i'fl - 2AeVi7'^/l - ^fViiFhF^fW - m2hF^hl'Fifl 
-20 flM'FhFf^^' - Ui'flf! - Uf^F^hj^Fif^ + 2flflF^-f^Fif2 - ^f^f^FF' f^-f' 
-3/2/iCVi7'/l^' ' 8fiF'fh'Ftfi - 12 fii^Fh^^Ffl - 2fUFfWF^fl - UfiM^F^^Ffl 
-2fUFhl^F^f^ ~ 20 fIfU'Fj'Ffi - Af^hiFF^fh^ - 32f2M'Fh^'fl + 2fihF'j*Fifl 
-8f2i^Ffl^^Ffl ~ 3.2flhF^fa^i^fl - ?,flhF^hi^F^f2) > , (52) 

where J has been defined in Eq. (pS)) . ft is interesting to notice that for these theories the reduction of the degrees of 
freedom does not occur as the Kasncr-typc background approaches the FLRW one, since det A does not either vanish 
nor diverge in the FLRW hmit. In fact, for the hmiting FLRW case we have 

Q F4-4Hf 

^^*^ = 4^nrTW^^^' ^''^ 

^ i2H\F + mHnF + mi) 

F-^{F + 2HF + 12m iY ' 

which arc both positive when F + 4iJ^ > as well as the odd-parity mode. 

A famous example is /(i?)-gravity. The result for this particular case is obtained by simply setting ^ = in the 
above equations. More explicitly, we have 

det A = 3^'(A + 2./2)^ 

[i^/1/2 + 2/i/2i^ + F/2/1 + 272/2(^/2 + /2F)]2 ' 

^ ^ ^F {Mfl + AVI + /2./i/i.f2 + 3 /|/27^/i + 3 /1/27VI + 3 /I7VI) ,5g^ 

[F/1/2+ 2/1/2^ + F/2A+ 27^/2(^/2 +/2F)]2 

In this case the sound speed is easy to evaluate, and one can verify that all modes propagate at the speed of light. 



C. Ghost crossing 

For the particular models discussed in this section, all modes become ghost when F + AH^ crosses 0. It might be 
interesting to study if this ghost crossing (i.e. the instant of time where any of the perturbation modes becomes a 
ghost) can appear without the background reaching a spacetime singularity. 

An example that shows ghost crossing can be easily constructed in the FLRW model. For dcfiniteness we consider 
models specified by f{R, G) = f{R + ^G)- In this case, the action can be written with the aid of one auxiliary field as 

S=^Jd^xV^[iR + \G)F- ViF)] . (57) 



167r 

Here ^ in the above calculation is to be identified with XF. Hence, the ghost crossing occurs when F + AXHF = 0. 
As long as the function V{F) can be chosen freely, the only non-trivial equation of motion is 

(1 + AXH^) F- H{1- 8XH + AXH^)F + 2HF = 0. (58) 

As a simple example, we assume power law expansion, a oct"'. Then, F is solved as 



F = wi 2Fi I 



. + 1 + Vn^+lOn+1 n + 1- \/n^ + lOn + 1 n + 1 t'^ 
4 ' 4 ' 2~'"4n2A 
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^«+3 TP M^ + 5 - Vn'^ + lOn + 1 n + 5 + Vn"^ + 10?i + 1 n + 5 \ 

2^^i(^ 4 , 4 '^'-^)' ^^^> 

with wi and W2 being constants. By choosing wi and W2 appropriately, we can make models which pass through 
F + AXHF = 0. For example, if we set n = 2 and W2 = 0, we have 

+ 64At2 + 256A2 , , 

F + AXHF = w, , (60) 

For negative A, this crosses when = — 16A/3 or — 16A. Since the values of R and Q are respectively given by 

6n(2n-l) 
i2 ' 

g= \, : (61) 

they are both regular at the ghost crossing pointijl. Here we discussed the simple FLRW case, but the ghost crossing 
without singularity generically occurs even if consider Kasner-type backgroun dfl 



V. CONCLUSIONS 



We studied the behavior of the perturbation for general modifications of gravity whose Lagrangian consists of a 
general function of R, the Ricci scalar and of G, the Gauss-Bonnet scalar, on Kasner-type backgrounds. We have 
shown that the existence of ghost is inevitable for generic models. However, the ghost mode can decouple from the 
physical spectrum in the FLRW limit. The kinetic term of this mode vanishes in this limit, which does not mean 
strong coupling here because the mass term does not vanish. Hence such cosmological models may have a chance to 
survive as an effective theory that describes only small deviation from the FLRW universe and/or the modifications 
of gravity tend to vanish at early times. We also found that the modified dispersion relation of the kind oj (x k'^ first 
discussed in Q is due to integrating out the ghost mode in the FLRW limit. 

Also, we presented special cases that avoid the presence of the ghost. In this case one degree of freedom is absent 
from the beginning since the action satisfies a special relation, and there is no reduction of degrees of freedom in the 
FLRW limit. If we require this condition to be satisfied for rather generic backgrounds, the form of the Lagrangian 
is restricted to £ = F{(j))R + ^{(j))G — y{<t>)- For such models, the ghost crossing, where the sign of the kinetic term 
flips for some modes, can happen. We have presented an example of ghost crossing within FLRW models without the 
background reaching a classical singularity. 
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